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摘要 : 望远镜 的 仪器 偏振 效应 主要 由 金属 eee. 是 当前 高 精度 天 文 偏 
振 测 量 主 要 面临 的 技术 问题 。 为 了 能 够 消除 仪器 偏振 对 偏振 测量 精度 的 影响 , 在 实验 室 对 金 
属 镀膜 材料 的 偏振 特性 进行 精确 而 全 面 测量 > 
种 偏振 特性 测量 方法 ， 该 方法 使 用 两 种 光电 调制 方式 对 保护 铝 膜 的 Mueller 矩阵 关键 元 素 进 
行 直接 测量 ， 并 基于 双 层 膜 的 结构 模型 ， 使 用 多 入 射 角 多 波长 下 的 测量 数据 在 不 依赖 色散 模 
型 以 及 参数 假设 的 情况 下 反 演 光学 参数 。 本 文 方法 的 反 演 结果 与 传统 方法 得 到 的 结果 保持 一 
致 ， 同 时 本 文 方法 能 够 对 偏振 特性 进行 更 好 的 表征 。 
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偏振 测量 是 天 理 物 理 研究 和 太阳 物理 研究 的 强大 观测 工具 ,利用 它 我 们 可 以 直接 获得 很 多 
天 文 目 标的 重要 物理 和 几何 信息 ， 而 这 些 信息 是 其 它 观测 手段 无 法 获得 的 。 例 如， 太阳 望 远 
镜 利 用 偏振 光谱 观测 诊断 太阳 大 和气 中 的 磁场 结构 m; 恒星 的 发 出 的 偏振 辐射 拥有 更 为 丰富 的 
物理 内 涵 叫 。 由 于 目标 的 低 偏振 特性 ， 天 文 偏振 测量 始终 都 是 一 种 光子 饥饿 的 观测 模式 。 因 
此 , 利用 大 口径 望远镜 收集 更 多 的 光子 ,已 提高 目标 偏振 信号 的 信 噪 比 就 成 了 天 文 偏振 测量 
的 重要 发 展 方向 。 望 远 镜 的 仪器 偏振 效应 ， 即 望远镜 光学 系统 对 入 射 偏振 信号 的 改变 ， 是 当 
前 高 精度 天 文 偏振 测量 主要 面临 的 技术 问题 。 对 于 大 口径 反射 望远镜 , 仪器 偏振 主要 是 由 金 
属 镜面 的 反射 偏振 效应 引起 。 根 据 菲 涅 尔 公 式 ， 金属 镜面 的 损耗 和 延迟 是 入 射 角 与 镀膜 材料 
折射 率 的 函数 。 从 微观 上 看 , 金属 膜 是 由 大 量 线 尺度 在 儿 十 纳米 的 金属 晶 粒 随机 堆 芭 而 成 的 ， 
其 光学 特性 是 镀 腊 材料， 镀膜 工艺 ， 使 用 环境 ， 以 及 使 用 时 间 的 函数 申 。 因 此 ， 为 了 能 够 准 


- 确 把 握 金 属 薄 膜 的 偏振 特性 ,进而 对 望远镜 进行 准确 的 偏振 光学 设计 或 建 模 ， 并 最 终 消 除 仪 
© 器 偏振 对 偏振 测量 精度 的 影响 , 我 们 需要 在 实验 室 对 金属 镀膜 材料 的 偏振 特性 进行 精确 而 全 
面 测量 研究 。 


在 大 口径 望远镜 中 常用 的 是 铝 膜 反射 镜 , 其 偏振 特性 已 被 广泛 研究 。 纯 铝 膜 镀 制 完成 不 久 
之 后 ， 会 在 其 表面 生成 天 然 的 氧化 铝 薄膜 &3。 氧 化 ea e er 但 也 
改变 了 铝 镜 的 偏振 特性 ， 且 在 表征 铝 镜 的 偏振 特性 的 时 候 ， 必 须 考 虑 氧化 膜 的 存在 加。 然而 
氧化 铝 薄膜 仍然 阻止 不 了 水 分 等 的 污染 , 因此 在 镀 完 铝 膜 后 会 继 乡 eee Nae 
硅 保 护 膜 B.7， 同 样 其 偏振 特性 会 有 很 大 的 变化 。 
传统 椭 偏 测量 方法 是 对 固定 入 射 角度 下 的 椭 偏 光谱 数据 (vw,A) 进 行 拟 合 ， 在 拟 合 过 程 中 会 
使 用 色散 模型 来 约束 各 个 波长 下 的 折射 率 E 0， 有 时 还 需要 使 用 已 知 的 文献 值 去 约束 参数 变 
化 。 因 此 ,传统 方法 的 反 演 的 参数 准确 度 很 大 程度 上 依赖 于 色散 模型 的 选取 ， 而 且 由 于 
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度 的 减少 与 参数 假设 等 的 限制 , 对 于 复杂 的 样品 需要 建立 十 分 复杂 的 结构 模型 才能 进行 较 好 
的 表征 。 

本 文 提 出 一 种 保护 铝 膜 偏振 特性 测量 方法 ， 该 方法 首先 对 保护 铝 膜 的 Mueller 矩阵 关键 元 
素 进行 直接 测量 ， 然 后 基于 双 层 膜 的 结构 模型 ， 使 用 测 得 的 Mueller 矩阵 数据 在 不 依赖 色散 
模型 以 及 参数 假设 的 情况 下 反 演 光学 参数 。 


1 测量 原理 与 方法 


天 文 目标 发 出 的 偏振 辐射 的 偏振 状态 可 以 用 Stokes 向 量 表示 , 即 Sn =O Q U V, 其 
中 工 表 示 入 射 光 总 强度 ，Q 为 垂直 与 水 平 偏振 方向 的 强度 参数 ，U 表示 土 45” 偏 振 方向 的 强 
度 参数 ,V 为 圆 偏 振 的 强度 参数 , 括号 外 的 上 角 标 T 表示 矩阵 转 置 。 经 过 镜面 反射 后 的 Stokes 
句 量 表示 Sout = MSn， 式 中 M ARARA ADERE, CRER RA ERU T: 


Mop Mọ»; My»; Myo) 
mo moso myso My 
M = 1>Q QQ U>Q VQ , #(1) 


Maou Maesu Myusy Myay 


May Mesy Myusv Myosy 


FP Ba ERE SCR ap ASR, FRAP da RA: I> XX = QUTY 表示 仪器 起 
偏 ，X1 > X241:X12 = QUTY 表示 仪器 串扰 ; X > X,X=Q, UV 表示 仪器 退 偏 。 

当 入 射 与 出 射 的 Stokes 向 量 按 如 图 1 所 示 的 偏振 坐标 系 ， 即 镜面 反射 变换 的 本 征 坐 标 系 
定义 时 ， 其 中 + Q 对 应 光 的 s 分 量 ， 金 属 镜面 反射 的 纪 勒 矩阵 可 以 简化 为 02]: 


ncoming light 


图 1 入 射 光 与 出 射 光 的 偏振 坐标 系 


Fig. 1 The polarization coordinate frames of incoming and outgoing light 


1 cos2W 0 0 
_ | cos2y 1 0 0 
M= | 0 0 sin2WoosA ins » #(2) 


0 0 —sin2WwsinA sin2ycosA 


FL pA ENG EBL FG WE Bl Sie JR ISIS EB E EIR E, A PEE BE 
Mog Myv mvav 三 项 表征 。 因 此 ， 在 随后 的 讨论 中 ，Mueller 矩阵 仅 指 代 以 上 三 个 非 零 阵 
元 。 


实验 数据 采集 是 使 用 商业 的 相位 调制 椭 偏 仪 (UVISEL-Plus, Horiba), 其 系统 结构 如 图 2 所 
示 . 从 光源 ( 握 灯 ) 发 出 的 光 经 过 一 个 起 偏 器 ( 洛 匈 棱镜 ) 后 斜 入 射 到 样品 表面 ,由 于 样品 的 偏振 
特性 ,反射 光 的 偏振 会 发 生 改变 ,之 后 经 过 光 弹 调制 器 PEM) 和 分 析 器 ( 洛 匈 棱镜 ) 的 调制 之 后 
被 光谱 仪 所 探测 ,最 后 得 到 三 个 测量 量 (10,16,10), 取决 于 基于 反射 面 的 起 偏 器 的 方位 角 P, 调 制 
器 的 方位 角 M, 分 析 器 的 方位 角 4 与 枯 偏 参数 (mA), 表 示 为 03]; 
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Ip = 1 + cos 2P cos 2M cos 2 (A — M) — {cos 2P + cos 2M cos 2(A — M)}cos 2 
I, = sin2(A— M)sin2Psin2 sin A (3) 
I, = sin 2(A — M) {sin 2M (cos 2W — cos 2P) + sin 2P cos 2M sin 2 cos A} 


图 2 ”相位 调制 椭 偏 仪 的 光学 结构 示意 图 


Fig.2 Optical configuration of the phase modulation ellipsometer 


BORK ATOR EIN EL BO RTE 1, | BOE RB BL, PFE 
偏 参 数 转化 为 样品 的 特性 。 TAT JECT She APPA Cini Ss BEE 17 Bz AR BIO IE SK PR EVE, 其 定义 
JHE BR Hi E[-1,1] HT, i ike, TARPS OGR, SRA. 
为 了 计算 的 准确 性 ， 我 们 可 以 直接 通过 对 到 ,大 数据 的 拟 合 来 研究 铝 镜 的 偏振 特性 。 由 式 (3) 
可 以 知道 通过 调整 偏振 器 ,调制 器 与 分 析 器 的 方位 角 , 我 们 可 以 得 到 不 同 的 测量 数据 .对 比 式 
CO) 可 知 , 铝 镜 的 Mueller 矩阵 的 几 个 非 实 数 元 素 可 以 由 两 种 调制 下 的 I, 表示 ,相应 测 得 数据 
如 下 : 


configuration 1:P = 45°,M = 0°,A = 45° 
Iz = sin2y sin A 
Ia = sin2WcosA 
configuration 2:P = 45°,M = 45°,A = 9 
1s2 = sin2y sin A 
I = cos 2Y 
此 时 , 我 们 就 可 以 直接 测量 到 Mueller MERE AI TERE, HH FAE BRATR TL 
可 以 得 到 反射 镜 更 多 的 信息 ， 则 可 以 把 测量 到 的 铝 镜 测量 数据 表示 成 Mueller 矩阵 形式 : 
1 Ie 0 0 
I2 1 0 0 
0 0 lo Isı 
0 0 =la Tet 


为 了 研究 铝 镜 偏振 特性 ， 我们 需要 获得 足够 多 的 信息 ， 鉴 于 仪器 测量 能 力 , RIKE TE 
反射 镜 样 品 (PF05-03-G01, THORLABS) 在 多 入 射 角度 (角度 范围 50-80”， 间 隔 2”) 下 的 
光谱 数据 (波长 范围 400 nm-1200 nm， 间 隔 25 nm)。 图 3 显示 的 是 测量 得 到 的 Mueller 矩阵 
关键 项 (mj >Q? Mysy» My CNT AY ee ES 代表 的 是 元 素 的 大 小 。 可 以 发 现 这 三 个 Mueller 
和 矩阵 元 素 的 光谱 随 着 入 射 角 的 变化 也 有 明显 变化 。 


0 H(A) 


#(5) 
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图 3 测量 得 到 的 Mueller 矩阵 关键 元 素 


Fig. 3 The key elements of measured Mueller matrix 


2 数据 分 析 


与 纯 铝 膜 相 比 , 加 了 氧化 硅 保 护 膜 的 铝 镜 的 偏振 特性 变 得 更 加 复杂 。 图 4 给 出 了 理性 情况 
下 带 有 保护 层 的 铝 膜 分 层 结构 ,其 中 mm 为 空气 的 折射 率 ,my 为 氧化 硅 的 折射 率 , ny 是 金属 铝 
的 复 折射 率 (定义 np = 7, 一 ik) ; 0m 为 在 空气 中 的 入 射 角 , 9p 是 氧化 层 内 的 折射 角 , 而 9 此 
时 是 复数 ,不 再 有 简单 的 折射 角 意义 ,但 仍然 服从 snell 定律 。 


1 
1 
1 
1 
ma 
1 
1 
i 


| 4 保护 铝 膜 的 理想 结构 模型 示意 图 


Fig.4 The ideal structural model of protected aluminum coating 


因此 符合 上 图 的 菲 涅 尔 反 射 系数 为 : 


(Mm — Np) cos d¢ +i (2-w) sin õp 
r= sv 7 #(6) 


Mm + Np) cos d¢ +i (2 +n) sind 
f 


其 中 6j = 至 nrdr cos 6/ 为 氧化 膜 的 相位 厚度 ,而 9 代表 有 效 折射 率 ， 其 下 标 ] = mn, fb 分 别 对 


应 于 空气 ,氧化 层 , 铝 膜 三 种 材料 ,而 s 与 p 偏振 分 量 的 有 效 折射 率 分 别 表示 为 : 
s-component: nj; = nj cos 6; 
ences nj = ny cos 6; mn) 
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根据 式 (2) 与 式 (6)， 并 结合 椭 偏 参数 的 定义 Ve = rj,/rs， 我们 可 以 为 图 4 中 的 结构 建立 偏 
振 模 型 : 


M = f(n,(A),K(A),n¢(A),dpstmAOm)> #(8) 


其 中 mh 、mr、r 与 nj 都 是 人 的 函数 。 
由 式 (9) 可 以 知道 偏振 特性 与 物理 参数 之 间 的 关系 ， 在 本 文中 假设 nm = 1, 对 于 特定 入 射 角 
与 特定 波长 下 的 Mueller ERER EDM g Myou Mysy DtFEMMAKAM T: 


es = fi(nmenpdy) 


Myoy = fo(n,K,npdp) #(9) 
Myou = f3 (nxnp,dy) 


若是 使 用 测量 到 的 Mueller 矩阵 元 素 ， 并 利用 上 述 关系 求解 光学 参数 ， 此 时 只 能 建立 三 个 

方程 。 本 文 的 方法 不 使 用 任何 文献 值 , 就 需要 求解 四 个 参数 , 方程 个 数 小 于 未 知 参数 的 个 数 ， 

无 法 求解 。 多 个 入 射 角 数据 在 单 层 铝 膜 的 时 候 并 不 能 得 到 更 多 的 信息 , 而 对 于 带 氧化 层 的 结 

构 ， 光 在 氧化 层 间 多 次 反射 之 后 ， 不 同 入 射 角 就 可 以 提供 额外 的 信息 [多 。 因 此 我 们 可 以 加 

~ 入 更 多 的 入 射 角度 下 的 测量 数据 以 增加 方程 ， 每 增加 一 个 入 射 角 度 ， 就 可 以 增加 三 个 方程 ， 

理论 上 只 需要 两 个 角度 就 可 以 实现 对 光学 参数 的 超 定 ,一 般 会 使 用 尽 可 能 多 的 方程 来 反 演 光 

学 参数 。 而 光学 参数 之 间 是 存在 相关 性 的 ， 例 如 与 x 是 正 相 关 na， 为 了 限制 其 相关 性 的 影 

响 ， 我 们 同时 拟 合 了 多 个 波长 下 的 数据 ， 每 个 波长 下 的 k 值 不 同 ， 而 dy 是 统一 的 。 此 时 在 拟 

合 光 谱 数 据 的 时 候 , 传统 方法 习惯 使 用 色散 模型 来 对 各 个 波长 下 的 折射 率 进行 约束 有 旦 只 使 用 

单个 入 射 角 的 数据 , 拟 合 效果 与 反 演 参数 准确 度 高 度 依赖 模型 的 选取 ,而 且 由 于 自由 度 的 减 

少 与 参数 假设 等 的 限制 , 复杂 的 样品 需要 建立 十 分 复杂 的 结构 模型 才能 进行 较 好 的 表征 。 本 

文 基于 图 4 所 示 双 层 膜 结构 模型 , 使 用 多 个 入 射 角 的 光谱 数据 , 同时 拟 合 多 个 波长 点 下 的 n、 

Ko nj 与 共同 的 dy， 且 不 使 用 色散 模型 ， 此 反 演 方法 自由 度 更 高 ， 能够 对 多 个 入 射 角 下 的 光 

谱 数 据 进行 更 好 的 表征 。 我 们 使 用 Levenberg-Marquardt 算法 L6 来 解决 这 个 非 线 性 最 小 二 乘 
优化 的 问题 : 


n 
l 2 
min) [Ce = msg) + Cer — mgu)? + (lsz 一 mysy)?|#(10) 
i=1 


其 中 n 为 铝 镜 Mueller 的 个 数 ; Mog? Mysu» My ,JJ 分 别 为 计算 得 到 的 关键 项 。 

入 射 角度 的 数据 越 多 能 提供 更 多 的 信息 , 反 演 未 知 参 数 的 精度 也 会 越 高 ， 因 此 我 们 会 用 上 
文 描述 的 方法 同时 拟 合 50-80” 入 射 角 下 的 光谱 数据 ， 并 使 用 传统 方法 反 演 结果 作为 对 比 。 
在 传统 方法 中 , 铝 膜 复 折 射 率 使 用 Drude-Lorentz 色散 模型 , 氧化 硅 使 用 Lorentz 色散 模型 然 
后 分 别 对 入 射 角 60° 与 70”〔 入 射 角 选 在 样品 的 腹 布 侍 斯 特 角 附近 ， 测 量 数据 灵敏 度 高 ) 
下 的 业 与 A 数据 进行 拟 合 。 

拟 合 的 效果 由 计算 的 值 与 原始 数据 的 残 差 与 原始 数据 的 比值 来 决定 , 记 为 相对 残 差 , 单位 
为 百分比 。 图 5(a) 所 示 是 两 种 方法 在 60° 与 70” 入 射 角 下 的 相对 残 差 ， 其 中 实 线 代 表 本 文 
方法 的 相对 残 差 ， 虚线 表示 传统 方法 的 相对 残 差 。 对 于 my， 0， Mysu» MysU» AS MC HE 出 的 
方法 都 比 传统 方法 的 拟 合 效果 好 ;尤其 是 mo 元 素 ， 传 统 方法 的 相对 残 差 到 了 50% 左 右 ， 
而 本 文 的 方法 的 相对 残 差 却 几 乎 为 零 。 本 文 方法 是 对 多 个 入 射 角 下 的 数据 进行 的 拟 合 ， 得 到 
的 拟 合 结果 是 符合 多 入 射 角 下 的 偏振 特性 ， 对 于 其 他 角度 的 相对 残 差 也 能 达到 相同 的 水 平 。 
此 外 , 两 种 方法 的 残 差 曲线 中 同时 出 现 相 对 残 差 突然 变 大 的 波长 点 ,这 是 因为 样品 在 这 个 波 
长 点 下 的 偏振 特性 特别 弱 ， 几 乎 为 零 ,〈 例 如 在 入 射 角 70” 下 ， 其 中 一 个 突变 点 725nm 处 的 
仪器 起 偏 特性 只 有 7e-4) 其 作为 分 母 求 得 的 相对 误差 都 会 变 得 特别 大 。 若 不 考虑 个 别 突变 点 ， 
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两 条 曲线 之 间 比 较 , 还 是 可 以 知道 本 文 方法 的 拟 合 效 果 更 好 。 表 1 中 显示 的 两 种 方法 反 演 的 
氧化 层 厚 度 ， 它 们 之 间 只 相差 0.5nm 左右 ， 两 种 方法 测 得 的 厚度 基本 一 致 。 


表 1 ”本 文 方法 与 传统 方法 反 演 的 dy 
Table 1 The retrieval d f of this work and traditional method 
This work Traditional(60° ) Traditional(70° ) 


d;/nm 190.5875 190.2589 190.0337 


两 种 方法 的 拟 合 的 光学 参数 与 文献 参考 值 如 图 5(b)(c)(d) 所 示 ， 其 中 实 线 是 由 多 个 波长 点 
的 结果 平滑 连接 而 成 ， 虚 线 是 铝 膜 复 折射 率 CORE McPeak) 与 氧化 硅 的 折射 率 (来自 
Gaol!) 的 文献 参考 曲线 。 对 于 mw 、r 与 nr, 其 图 中 三 条 曲线 的 变化 趋势 基本 保持 一 致 ， 只 是 
本 文 方法 与 传统 方法 在 数值 上 有 一 点 差别 , 一 是 因为 传统 方法 使 用 色散 模型 并 不 能 完全 表征 
反射 镜 的 偏振 特性 ， 其 中 忽略 了 很 多 影响 偏振 特性 的 因素 ， 二 是 因为 本 文 方法 考虑 了 mi,o 
的 变化 特性 。 尤 其 对 于 图 5(b)， 本 文 工 作 与 虚线 即 参考 曲线 在 大 于 1000nm 波段 的 差异 明显 
增 大 ， 参 考 曲 线 也 是 由 传统 方法 求 得 ， 也 是 在 70” 入射 角 下 对 他 的 铝 膜 样品 进行 的 测量 。 
因此 可 以 发 现 参考 曲线 与 传统 方法 在 70” 下 测量 本 文 的 样品 的 十 分 相近 ,但 在 大 于 1000 nm 
的 波段 下 还 是 有 所 差异 。 对 于 传统 方法 在 入 射 角 60” 下 的 结果 ， 同 样 在 大 于 1000 nm 的 波 
段 与 参考 曲线 有 所 差异 。 因 此 本 文 工 作 与 参考 曲线 的 差异 可 以 认为 是 样品 的 差异 导致 。 而 本 
文 工 作 与 传统 方法 曲线 之 间 的 差异 可 以 认为 是 不 考虑 样品 各 向 异性 的 前 提 下 , 均衡 多 个 入 射 
角 下 偏振 特性 的 结果 。 本 文 方法 的 曲线 都 不 够 平滑 , 是 因为 本 文 方法 没有 使 用 色散 模型 等 方 
式 对 光学 参数 随 波 长 的 变化 趋势 进行 约束 , 只 使 用 共同 的 氧化 膜 厚 度 去 建立 各 个 波长 点 之 间 
的 联系 ， 得 到 的 结果 反而 能 够 真实 反应 当前 波长 点 下 的 偏振 特性 。 
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图 5 (a) 两 种 方法 的 在 60” 与 70” 入 射 角 下 对 于 测量 数据 的 相对 残 差 ; 
(Cb)(c)(d) 分 别 为 两 种 方法 反 演 的 nm 、r* 与 mr 随 波 长 变化 的 对 比 图 


Fig. 5 (a) The relative residuals for measured data by this work and traditional method at 60° and 70° ; 


(b)(c)(d) is the comparison of results for the dependence of wavelength of n,, K, np, respectively 


本 文 首先 介绍 了 对 保护 铝 膜 的 Mueller 矩阵 关键 元 素 进行 直接 测量 的 原理 与 方法 ， 然 后 介 
绍 了 在 不 依赖 色散 模型 以 及 参数 假设 的 情况 下 反 演 光学 参数 的 方法 与 结果 。 实 验 结果 表明 ， 
该 方法 反 演 得 到 的 光学 参数 与 传统 方法 的 结果 一 致 。 相 比 于 传统 方法 , 该 方法 获取 了 更 多 的 
偏振 特性 信息 ， 且 自由 度 更 高 ， 能 够 对 偏振 特性 (尤其 是 仪器 起 偏 特性 mj g) 进行 更 好 的 
表征 。 本文 方法 的 结果 曲线 都 不 够 平滑 ,因为 本 文 方法 没有 使 用 色散 模型 等 方式 对 光学 参数 
随 波 长 的 变化 趋势 进行 约束 , 得 到 的 结果 反而 能 够 真实 反应 当前 波长 点 下 的 偏振 特性 。 本 文 
方法 是 对 多 个 入 射 角 下 的 数据 进行 的 拟 合 ， 得 到 的 反 演 结果 是 符合 多 入 射 角 下 的 偏振 特性 ， 
而 且 基 于 简单 的 双 层 膜 结构 ， 更 有 利于 未 来 对 未 知 入 射 角 下 的 偏振 特性 进行 预测 与 验证 
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Abstract: The instrument polarization effect of the telescope, which is mainly caused by the 
reflection polarization effect of metal coatings on mirrors, is the principal technical problem 
facing the current high-precision astronomical polarimetry. It is essential to perform accurate and 


comprehensive measurements and research on the polarization properties of metal coating 


materials in the laboratory in order to eliminate the influence of instrument polarization on the 
accuracy of polarimetry. This paper proposes a polarization properties measurement method for 
the protected aluminum coating. This method can directly measure the key elements of the 
Mueller matrix of the protected aluminum coating by two photoelectric modulations at a 
multiple-angle-of-incidence spectrum. This method also retrieves the optical parameters based on 
the double-layer model without dispersion model and parameter assumptions. The retrieval results 
of the method in this paper are consistent with those obtained by traditional methods, and the 
method in this paper can fully characterize the polarization properties. 

Key words: Polarimetry; Instrument polarization; Polarization properties of protected aluminum 
coating; Mueller matrix; Retrieve of optical parameters 


